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. The fluctuations of the Higgs and pseudo Nambu-Goldstone fields in the 2SC phase with mis- 

' matched pairing are described in the nonhnear reahzation framework of the gauged Nambu-Jona- 

Lasinio modeL In the gapless 2SC phase, not only Nambu-Goldstone currents can be spontaneously 
generated, but the Higgs field also exhibits instablity. The Nambu-Goldstone currents generation 
indicates the formation of the single plane wave LOFF state and breaks rotation symmetry, while the 
Higgs instability favors spatial inhomogeneity and breaks translation invariance. In this paper, we 
focus on the Higgs instability which has not drawn much attention yet. The Higgs instability cannot 
be removed without a long range force, thus it persists in the gapless superfluidity and induces phase 
separation. In the case of g2SC state, the Higgs instability can only be partially removed by the 
^ , electric Coulomb energy. However, it is not excluded that the Higgs instability might be completely 

OO ■ removed in the charge neutral gCFL phase by the color Coulomb energy. 
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I. INTRODUCTION 



Color superconductivity (CSC) 0,0] of QCD at high baryon density has been an active area of research since 1998. 
" It involves both relativistic field theory and statistical mechanics, and many novel properties of this unusual phase 
of QCD have been revealed. For reviews on recent progress see, for example, Ref. The color superconducting 
^ ' phase is expected to reside in the central region of compact stars, where the charge neutrality condition as well as /3 
• ^ , equilibrium are essential |3| ■ 

' The research of charge neutral cold dense quark matter has driven the color superconductivity theory into a new 
territory beyond standard BCS theory. The requirement of charge neutrality condition induces a substantial mismatch 
between the Fermi surfaces of the pairing quarks, which reduces the available phase space for Cooper pairing. The 
results based on conventional BCS framework brings us puzzles. For example, in the standard BCS framework, it 
was found that at moderate mismatch, homogeneous gapless superconducting phases can be formed. However, 
gapless superconducting phases exhibit anti-Meissner screening effect, i.e., chromomagnetic instability 2^ 8], which is 
in contradict with the Meissner effect in standard BCS superconductivity. 

Superfluidity or superconductivity with mismatched Fermi momenta also appears in other systems, e.g., electronic 
superconductor in a strong external magnetic field |lCl lllj , asymmetric nuclear matter |12| , and in imbalanced cold 
atomic systems It has been an unsolved old problem how a BCS superconductivity is destroyed as the 

mismatch is increased and what exotic phases may occur before entering the normal phase. 

It was proposed in the 1960s that, applying a strong external magnetic field in an electronic superconductor, 
the competition between the pair breaking and pair condensation would induce an unconventional superconducting 
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phase, i.e., Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) state In the LOFF state, the Cooper pair carries a net 

momentum g, which is different from the zero net momentum BCS Cooper pair. The order parameter of the FF state 
is characterized by a single plane wave: A e^^'^, and the LO state is a superposition of two plane waves or "striped 
phase" : A cosg • x, where A is the amplitude of the gap. In general, the structure of the order parameter can be more 
complicated. The LOFF state has still not yet been confirmed experimentally in electronic superconduting systems 
[la, llal due to the difficulties from orbital effects and impurities, and it still remains being persued after more than 
40 years. 

Imbalanced cold atom systems offer another intriguing experimental possibility to understand how Cooper pairing 
is destroyed. By making the populations of atoms in two different hyperfine states unequal, one controls the mismatch 
between their Fermi surfaces. Due to the absence of both the orbital effects and impurities, it seems very promising to 
search for the LOFF state in imbalanced cold atom systems. It has stirred a lot of theoretical interests of discussing the 
formation of the LOFF state in mismatched cold atom systems 0|. However, recent experiments |0 in imbalanced 
cold atom systems did not show evidence of the formation of the LOFF state, rather indicated a non-uniform state 
of phase separation state. 

The momentum-carrying Cooper pairing state was extended to color superconductor and was further discussed 
in Refs. |20|. However, we started to appreciate the importance of this unconventional superconducting state only 
after relating it to chromomagnetic instabilities of gapless CSC ,21]. Recent studies have demonstrated that the single 
plane wave FF state is unavoidably induced by the instability of the local phase fluctuation of the superconducting 
order parameter, i.e., the Nambu-Goldstone boson fields 22, 23, 24, 25, 26] . 

Obviously, there should be something missing for our understanding of the pairing breaking state if we failed to 
observe the (LO)FF state in the mismatch regime where the magnetic or superfluid density instability [13. l2^ develops. 
In Refs. [2a ] and 29], a new type of instability, Higgs instability, was discovered. The Higgs instability, being an 
instability with respect to an inhomogeneous fiuctuation of the gap magnitude, may favor non-uniform mixed state 
or LO state more than the single-plane wave FF state or Nambu-Goldstone current state. 

The non-uniform mixed phase has been extensively discussed as one possible candidate of the ground state of the 
asymmetric fermion pairing system 's^ by comparing the free energy, however, there was no much effort of trying to 
understand the origin of forming a mixed phase. We try the first step toward this direction by exploring the Higgs 
instability of a homogeneous 2SC. The Higgs instability is a straightforward extension of the Sarma instability but has 
not received much attention until recently '2^ [2^ . The Higgs instability refers to the instability of a homogeneous 
CSC state against an inhomogeneous fluctuation of the magnitude of the gap parameter, which is not protected by 
globally imposed constraints, its salient features have been discussed in the short letter [2^ . 

In this paper, we shall mainly display the details behind the arguments made in [29l ]. In Sec. ^ we describe the 
nonlinear realization framework of the gauged SU(2) Nambu~Jona-Lasinio(NJL) model in /3-equilibrium. In Sec. IIIII 
we investigate the phase fiuctuation as well as the amplitude fiuctuation, and show that with the increase of mismatch, 
the Nambu-Goldstone currents are spontaneously generated, and the Higgs field also becomes unstable. In Sec. IIVI we 
discuss the general features of the Higgs instability. The Higgs instability will naturally induce spatial inhomogeneity, 
however, whether phase separation can be developed in the system has to be determine by its cost of Coulomb energy, 
which is discussed in Sec. We give summary in Sec. IVII 



II. THE GAUGED SU(2) NAMBU-JONA-LASINIO MODEL 
A. The Lagrangian 



To analyze the Higgs instability of the g2SC phase, we start from the gauged Nambu-Jona-Lasinio (gNJL) model, 
the Lagrangian density has the form of 



C 



j(lp + At") q + Gs [{qqf + [qil^'Tq) 



G 



D 



q^i-i^T2e''q 



qi-fTitPq^ 



(1) 



with D,, 



9. 



igA'^T"'. Here A" are gluon fields, — A"/2 are the generators of SU{3)c gauge group with 
a = 1, • • • , 8. In the gNJL model, the gauge fields are external fields and do not contribute to the dynamics of the 
system. The property of the color superconducting phase characterized by the diquark gap parameter is determined 
by the nonperturbative gluon fields, which has been simply replaced by the four-fermion interaction in the NJL model. 
Gs and Go are the quark-antiquark coupling constant and the diquark coupling constant, respectively, g*-^ = Cif" , 
(f — q^C are charge-conjugate spinors, where C = 17^7^ is the charge conjugation matrix (the superscript T denotes 
the transposition operation). The quark field q = qia with i — u,d and a = r, g, 6 is a fiavor doublet and a color triplet, 
as well as a four-component Dirac spinor. r = (r^, r^, r^) are Pauli matrices in the fiavor space, with antisymmetric, 
and (e^)"^ is antisymmetric tensor in color space with p = r,g,b. fi is the matrix of chemical potentials in the color 



3 



and flavor space. In /3-equilibrium, the matrix of chemical potentials in the color-flavor space /t is given in terms of 
the quark chemical potential fi, the chemical potential for the electrical charge fj.^, and the color chemical potential 

A'S'' = - ^^eQ^J)S''^ + ■^^^8S^, (Tg)"^. (2) 

In this paper, we shall focus on the color superconducting phase, and neglect the influence from the chiral con- 
densates by assuming < qq and < qj^rq >= 0. The is the case when the baryon density is sufficiently high. 
Introducing the anti-triplet composite diquark fields 

AP{x) ^2iGDiq^t-f^T2ePq), A*P{x) ^ ~2iG niqil'' T2e'' q"^ ) , (3) 

one obtains the bosonized version of the model for the 2-flavor superconducting phase, 

C2SC = q{^P + At")? - 7^A*''A'' 

.[fi^'^T2ePq]-—[qil'^T2ePq^]. (4) 



In the Nambu-Gor'kov space, 



the Lagrangian Q simplifies to 



* = ( „c ) , (5) 



C2sc = -^^ + ^^[S]-^^. (6) 



The inverse of the quark propagator in momentum space is defined as 



1^ A+ [G, 



(7) 



with the off-diagonal elements 

A- = -^iT2eP^zA\ A+ = irse^TsA*", (8) 
and the free quark propagators (P) taking the form of 

[G^(^r'=7°(Po±A)-7-P, (9) 

with the 4-momenta denoted by capital letters, e.g., P = {po,p). In terms of the Pauh matrices with respect to NG 
indexes, 

the inverse propagator |7J can be cast into a compact form 

[S{P)]-^ = 7"(Po + P3A) -1-V- P+A- 4- p_ A+ (11) 

with p± = {pi ± ip2)/2. 



B. Nonlinear realization framework of the gNJL model 



The superconducting state is characterized by the order parameter A(x), which is a complex scalar field and 
has the form of A(x) = |A(a;)|e*''''-^'', with |A| the amplitude and ip the phase of the gap order parameter. For a 
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homogeneous condensate, A{x) is a spatial constant. The fluctuations of the phase give rise to the pseudo Nambu- 
Goldstone boson(s), while that of the amplitude to the Higgs field, following the terminology of the electroweak theory. 
Stimulated by the role of the phase fluctuation in the unconventional superconducting phase in condensed matter, 
we follow Ref. '33^ to formulate the 2SC phase in the nonlinear realization framework in order to take into account 
naturally the contribution from the phase fluctuation or pseudo Nambu-Goldstone current (s). 

In the 2SC phase, the color symmetry G = SU{3)c breaks to H = SU{2)c- The generators of the residual SU{2)c 
symmetry H are {S" = T"} with a = 1,2,3 and the broken generators {X'' = r^+S} with 6 = I,-- - ,5. (More 
precisely, the last broken generator is a combination of Tg smd the generator 1 of the global U{1) symmetry of baryon 
number conservation, B = {1 + ^/3Ts)/3 of generators of the global U{1)b and local SU{3)c symmetry. ) 

The coset space G/H is parameterized by the group elements 



V(x) = exp 



(12) 



here operator V is unitary, and V^^ = and (pa{a = 4, • • • ,7) and ips are five Nambu-Goldstone diquarks, and we 
have not consider the topologically nontrivial case and therefore V(a;) can be expanded uniformly according to the 
powers of (/s's. In fact, V{x) is alway topologically trivial for a configuration of (p^s that has a finite energy because of 
the trivial homotopy group Tr2{SU{3)/ SU{2)) pj . 

Introducing a new quark field x, which is connected with the original quark field q in Eq. (0J through a nonlinear 
transformation, 

9 = Vx , q^xV^ , (13) 

and the charge-conjugate fields transform as 

qc = V* xc , qc = Xc . (14) 

The advantage of transforming the quark fields is that this preserves the simple structure of the terms coupling the 
quark fields to the diquark sources, 

qcA+q = xc^^X, qA-qc = x^~XC- (15) 
In the Nambu-Gor'kov space of the new spinors 

^={xc) ' ^^(^'Xc), (16) 

the nonlinear realization of the original Lagrangian density Eq. (|SJ) takes the form of 

$+$- 1 _ 

^2SC^-^^ +^XS-,'X, (17) 

with 

Here the explicit form of the free propagator for the new quark field is 

[GLir' = ^^ + /i7o+7p^^ (19) 

and 

[GomV' = ip^-filo+Jf.VS. (20) 

Comparing with the free propagator in the original Lagrangian density, the free propagator in the non-linear realization 
framework naturally takes into account the contribution from the Nambu-Goldstone currents or phase fluctuations, 
i.e., 

= vt {id'') V, 

= {id'') V*, (21) 



] • (18) 
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which is the NcNj x iVc A^/-dimensional Maurer-Cartan one-form introduced in Ref. [s^. The Hnear order of the 
Nambu-Goldstone currents and has the exphcit form of 

8 

^ -^(aVa)r, , (22) 

a=4 

VS 0, Y^id^^a) {Taf ■ (23) 

The advantage of the non-Hnear reahzation framework Eq. (|17f) is that it can naturaUy take into account the 
contribution from the phase fluctuations or Nambu-Goldstone currents. The task left is to find the correct ground 
state by exploring the stability against the fluctuations of the magnitude and the phases of the order parameters. The 
free energy ^{Vfj,, $, /i, /is, /ie) can be evaluated directly and it takes the form of 

^lnl{v^,^,^i.^is.^ie) = -\TY J (0?Trin([5,4P)]-i) + (24) 



III. THE FLUCTUATION OF THE HIGGS AND NAMBU-GOLDSTONE FIELDS 
A. Expansion the free-energy in terms of the BCS Cooper-pairing state 

In the 2SC phase, the color symmetry SU{3)c is spontaneously broken to SU{2)c and diquark field obtains a 
nonzero expectation value. Without loss of generality, one can always assume that diquark condenses in the anti-blue 
direction, i.e., only red and green quarks participate the Cooper pairing, while blue quarks remains as free particles. 
The ground state of the 2SC phase is characterized by (A^) = A, and (A^) = 0, (A^) = 0. 

Considering the fluctuation of the order parameter, the diquark condensate can be parameterized as 




exp 



YVa{x)Ta 







A + H{x) 



(25) 



where $''(x) = (0, 0, A -I- H{x)) is the diquark field in the nonlinear realization framework, tpa,^?, are Nambu- 
Goldstone bosons, and H is the Higgs field. 

Expanding the diquark field around the ground state: = (0, 0, A), the free-energy of the system takes the 
following expression as 



NG 



H- 



(26) 



There are three contributions to the free-energy, the mean-field approximation free-energy part ^Im which has the 
form of 



T 



M 



E 



2 (27r)3 

the free-energy from the Higgs field has the form of 



^'P^Trln([5M(P)]-^) ' 



(f'k 

(2^ 



H*{K)nH{K)H{K) 



(27) 



(28) 



with 



T 



2G, 



E 



(2^ 



Tr 



SMiP + K) 







Sm{P) 



iT2e^j5 
-«T2e^75 



(29) 



and the free-energy from the Nambu-Goldstone currents ^ng has the form of 



n 



T 



NG 



4 EEy (27r)3 (27r)3 



Tr 



\Sm{P- 



i^^{~K)^^ 







(30) 
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with 



'>^iK) = gA^:{K)Ta-V^{K) , (31) 



col^iK) = -gA';:{K) Tj + V^iK) . (32) 



where the inverse propagator Sj^j^ takes the form of 

V -^T-se^TsA [Go {P)\ J 
= 7°(Po + P3A)-7-P + Ap2T2e^75- (33) 
The quasi-quark propagator at mean-field approximation has the form of 

Sm = ) . (34) 

and its expUcit expression of the Nambu-Gorkov components of Sm has been derived in Ref . • The weak couphng 
approximation of Sm can be found in the appendix A of this paper. 

The Matsubara self-energy functions Hh (K) of the Higgs field and that of the Goldstone fields (obtained after the 
sum over po and integral over p in Eq. H3()|l) can be continuated to real frequency following the standard procedure. 
Because of the gapless excitations, the values of these functions at zero frequency and zero mometum depends the order 
of the limit. In this work, we shall restrict our attention to static fluctuations only which amounts to replace H{K), 
A{K) and Lp{K) of Eqs. (EHJ and ^ by VTH{k)Sko,o, VTA{k)Skofi and VTip{k)Sko.o- Thus the long wavelength 
limit of the self-energy functions discussed below corresponds to the limit lim^^^^ limfeQ^o- 



B. Free-energy at mean-field approximation and Sarma instability 

In the mean-field approximation, the free-energy for u,d quarks in /3-equilibrium takes the form 

Here we did not take into account the contribution of free electrons. The sum over A runs over all (6 quark and 6 
antiquark) quasi-particles. 

At zero temperature, the mean-field free-energy has the expression of 

A2 A4 ,,4 ,,4 



M 



AGd 27r2 127r2 127r2 
'•^ p'dp 



^ (Vb + + A2 + ^{p - nY + A^) 

29 iSfi - A) r ^ (5^ _ V(p - fir + A2) . (36) 



As we already knew that, with the increase of mismatch, the ground state will be in the gapless 2SC phase when 
A < 6fi, the thermodynamical potential of which is given by 



7r2 V" Ao 

dynamic potential 
mismatch, (5/i = 0. The solution to the gap equation reads 



- ^ ^ In St^+ySj^'-^' _ S^^VSJ^^^ + (37) 



where fl^^ is the normal phase thermodynamic potential. Aq the solution to the gap equation in the absence of 



A - v/Ao(2<5/^- Ao). (38) 

The gapless phase is in principle a metastable Sarma state |ir)l |. i.e., the free-energy is a local maximum with respect 
to the gap parameter A. We have 

V aA2 J -,,5,- 7r2 V ^sn-^-A-^y ^ ' 

The weak coupling approximation is employed in deriving Eqs. H37|) and H39|) from Eq. (|36|) . which assumes that Ao, 
A and Sfi are much smaller than fi and A — /i. The same approximation will be applied throughoutthe paper. 
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C. Nambu-Goldstone currents generation and the LOFF state 



The quadratic action of the Goldstone modes in the long wavelength limit can be written down with the aid of the 
Meissner masses evaluated in Ref. 7j. We find that 



1 f ^ ^ ^1 ^1 
^ - / (fry^ "mliA."- V(p'')(A'' Vip'') + higher orders . 



where mi — m2 = — 0, 



2_ .gV [^'-2,5^*2 



37r^ 



2A2 



0{6n - A) 



and 



m 



9 -2 
2 _ 3 M 
8 



97r2 



6^9(3^ - A) 



(40) 



(41) 



(42) 



It was found that at zero temperature, with the increase of mismatch, for five gluons with a = 4,5, 6, 7, 8 correspond- 
ing to broken generator of SU{3)c, their Meissner screening mass squares become negative ?!. This indicates the 
development of the condensation of 



8 1 

^ < A'' - - V(p° >^ 0. 



(43) 



It can be interpreted as the spontaneous generation of Nambu-Goldstone currents J2a < ^V'" >7^ or gluon 
condensation X]a=4 ^ !3- ^^"^ interpreted as a colored-LOFF state with the plane- wave order 

parameter 



A(x) = Ae^^"=-'^'^°-^. 



(44) 



D. Higgs instability 

In the previous two subsections, we discussed the two known instabilities induced by mismatch, i.e., the Sarma 
instability and chromomagnetic instability, respectively. In this subsection, we focus on another instability, which is 
related to the Higgs field, and we call this instability "Higgs instability". 

The free-energy from the Higgs field can be evaluated and takes the form of 

= 2 ^i?*(fc)H«(fc)F(fc). (45) 



0-- 



FIG. 1: The Feynman diagrams for the Higgs polarization function IIh. 
Evaluating the one-loop quark-quark bubble H/f (fc) shown in Fig. ^ we obtain that 

i2 n2 



H^(A:) = ^Iik\S^l) + ^k'j{k\5y), (46) 
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where the functions I{k\S^) and J{k\5ii) are given by 



J(fc|(5/z) ^ A'rV , / dx- (48) 



with Re-\/ + i(5/i)2 + > 0. Here, the summation over n is the frequency summation at finite temperature field 
theory, with a;„ = (2n + l)7rT. At T=0, both functions can be simphfied as shown in the appendices A and B. Their 
power series expansion yield 

IlHik)^AH+BHk'' (49) 

with 



Bh 



2/22 
97r2A2 



(51) 



It follows from Eqs. H5()(l and (|51(l that the Higgs field becomes unstable in the gapless phase when (5/i > A, Ah in 
the gapless phase is shown in Fig. 01 by the red dash-dotted line. The Higgs instability was also considered in Ref. 
[2^ , where it was called as " amplitude instability" . It has to be pointed out that we got different expressions for the 
coefficients of the gradient term. For k » A, we have 



nff(fc) 



27r2 



k 
A 



(2 In ^ - 2 - In 



5^1 - y^Sj? 



6fi + \/6fJ^ 



(52) 



Therefore the Higgs instability disappears for sufficiently large momentum. The form factor Hij(fc) for arbitrary 
momentum is plotted using red dashed line in Fig. ^for a typical value of the mismatch parameter A/S^ = 1/2. We 
notice that the Higgs instability becomes stronger for nonzero momentum. 

Near the transition temperature, A << T, we may drop A's inside the square roots and in the denominators on 
RHS of Eqs. (I47|l and H48(l . The Higgs self-energy function (|46ll goes to the expression derived from the Ginzburg- 
Landau theory |2J|- The Ginzburg-Landau free energy for g2SC can be easily obtained from that of an electronic 
superconductor by replacing the Fermi momentum kp with fl, the Fermi velocity with one and multiplying the whole 
expression by four (four pairing configurations per momentum for g2SC versus one for the electronic SC). 

We would like to make two remarks before concluding this section: 1) The combination of A — |V((5 of Eq. 14U|) is 
by no means generic for a nonrelativistic system. The general expression of ^Ing beyond the long wavelength limit 
takes the form 



^NG — 7: 
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d^k 



(27r)3 ^ 



H^(fc)A'r(A?) • A"' (fc) + n\Uk){A?r + -kip''*{k)) ■ {A?,{k) - -kip^ik)) + higher orders 

"5 "5 



(53) 

in momentum representation, where A||(A^) is parallel(perpendicular) to k. Except that H;]"(0) = Ha(0), they may 

be quite different for fc ^ 0, so are the instability domains with respect to k. 2) At T — 0, the dependence of the 
thermodynamic function on the diquark field is highly nonlinear and the gradient term is not Ginzburg-Landua like, 
i.e., VA''* • VA''. There are no simple relations between the coefficient Bh and the Meissncr mass squares. This may 
explain the difference between our Bh and that of . 



IV. A GENERAL VIEW OF THE HIGGS INSTABILITY 



In this section, we shall provide a general view of the Higgs instability in relation to the Sarma instability as well as 
the Coulomb energy in relation to the charge neutrality constraint. Consider a system described by the Hamiltonian 
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T-L with a spontaneous symmetry breaking. The thermodynamic potential is given by 

f] = -rinTrexp ' ' (54) 

where Q's are m conserved charges with /i's the corresponding chemical potentials. The trace extends to a subset of 
the Hilbert space specified by n order parameters Aq with a = 1, rt, which are invariant under the symmetry group 
of Ti and are homogeneous. In the mean field approach, one approximates Ti with TLmf that depends on the order 
parameters explicitly and leaves the trace unrestricted. In the superfluid of imbalanced atoms or in g2SC, we have 
m = 2 and n = 1. Qi and Q2 stand for the number of atoms of each species in the former or the electric charge and 
the 8th color charge in the latter. Ai corresponds to the magnitude of the gap parameter in both cases. In particular, 
Ai = A of the last two sections. For gCFL, to = n = 3 with Qi, Q2 and the electric, 3rd color and 8th color 
charges while Ai, A2 and A3 correspond to the magnitude of the eigenvalues of the 3x3 diquark condensate matrix 
<j)C ^ gccic2g^^^^_^ ^ 9/^9/2 ^- "^^^ equihbrium conditions read 



and the expectation value of the Q-quanta density reads 



= ^<Q.>--(|^) (56) 

V va^j / Ai,...,A„ 



where V is the volume of the system and the thermal average < Qj > is defined as 



liQj exp I ^ 

< >= — ^n^^'^u,,o.^' ■ (57) 



IT exp I 

The Sarma instability corresponds to the negative eigenvalues of the n x n stability matrix 



Under the constraints of fixed densities rij 's, the relevent thermodynamic quantity to be minimized is the Helmholtz 
free energy which is obtained by a Legendre transformation, 

J=' ^n + ^HjUj. (59) 
j 

While the equilibrium conditions 

(S- ..^"^ <™' 

are equivalent to Eq. H55|l . the stability matrix becomes 



where the m x m matrix 



(^) <«2) 

\0^i0^j / Ai,...,A„ 

is positive. To see this point, we carry out the second order derivative and obtain that 

My =< Q,Qj > - <Qi >< Qj >=< {Q,- < Qi >){Qj- < Qj >) > . (63) 
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The positivity follows from the observation that Mz > for an arbirary m x 1 complex column matrix z. Therefore, 
the stability matrix (|61|) could be positive even Eq. (|58|l is not so. This is precisely what the constraint of fixed particle 
numbers in case of BP or the constraint of the charge neutrality accomplished. Because of the relations ()5fi|l and 

A^. = i^) > (64) 

we have 



M,,= ('|^) (65) 



■3 / Au 



and 



Therefore the positivity of Mij also implies that of the matrix 



iM-%^(P-) -(^] ■ (66) 



(67) 



The inclusion of inhomogeneous variations renders the stability matrix infinitely dimensional with additional matrix 
elements 



) = ( ^ (68) 



^(90*(A:)906(fc) ^"i,---,"™ ^d4>l{k)d4>b{k) ' Mi:---./' 

at fc 7^ 0, without the second term of Eq. (|61|l . where we have parametrized a general variation of the order parameters 
according to 

b^a + \Y.^i>a{ky^-' (69) 

with V the volume of the system. If the matrix elements Eq. (|68|) are not singular at fc = in the thermodynamic 
limit, V oo, 

1. I d'^^ \ I d'^n \ 

hm lim — = 70) 

fe^o^^-*°° ^a(/)*(/c)a(/)6(A:)^A'i,-,M„. vaAaaAfc/Mi.-.M™ 

The Sarma instability will show up with respect to the variations of a sufficiently low but nonzero momentumn. 

The Higgs instability can also be understood intuitively. Let us devide the whole system being considered into J\f 
subsystems of equal volumes and equal chemical potentials. We label them by the integer p. As long as the size of 
each subsystem is much larger than all microscopic length scales, the energy of the interfaces can be ignored. At 
equilibrium, the thermodynamic potential density of the master system is 

f^*°'- = -;^^r!(P)=r!M (71) 

p 

where J^m is given by (|37|l with A by H38|) . Now we consider an inhomogeneous variation of the gap parameters, 
SA^^\ 6A^n \ such that 1) each of the 7i-component vector, (A^^-*, A^^^) proportional to the same eigenvector of 
the matrix 



(72) 



\dAadAb 

with the negative eigenvalue A < 0, 2) J2p ^^a"^ — 0. Then 

^ i_ y |^<5Aip) = i- y 1^ y sM^^ = o (73) 

^ Af ^ dAa M ^ dAa ^ " ^ ' 

a,p 
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but 

Sn^ot. ^ 1 y 9^ y ^^(p)^^(rt ^ 1 A y 5 A(f )5 A(f ) < 0. (74) 

a,h p a,p 

Therefore the Sarma instabihty develops through different variations of the order parameters of each susystem while 
maintaining the constraint conditions of the master system. 



V. COULOMB ENERGY 



Unless a competing mechanism that results in a positive contribution to the inhomogeneous blocks of the stability 
matrix. The Higgs instability will prevent gapless superfluidity/ superconductivity from being implemented in nature. 
In the system of imbalanced neutral atoms, such a mechanism is not likely to exist and this contributes to the reason 
why the BP state has never been observed there. For the quark matter being considered, however, the positive 
Coulomb energy induced by the Higgs field of electrically charged diquark pairs has to be examined. 

The second order derivative of the Hemlholtz free energy, Eq. (|61|l for g2SC reads: 



KdA^Jn, ~ Ua2 



and the stability implemented by the charge neutrality implies that 

While the global charge neutrality is maintained, an inhomogeneous (f) will induce local charge distribution. The 
corresponding Coulomb energy 



_ 1 ^ spjkrspjk) 



should be added to the RHS of the stability matrix element H68|) where Sp{k) is the Fourier component of the ^-induced 
charge density and m^i is the Coulomb polarization function ( Debye mass at fc = ). We have 

^lik) = -^J2^noQSiP + KhoQSiP) (78) 
^ p 

and 

5p{k) = K{k)H{k), (79) 

where 

= ^ E tr7oQ5(P + K)j5e^p2SiK) (80) 
^ p 

with K = (0, k) and Q the electric charge operator. We have 

Q^P3{a + bT3) (81) 

with a = 1/6 and b = 1/2 for the quark matter consisting of u and d flavors. The Eq. H77|l corresponds to the 
diagram a of Fig. [21 The diagram b in Fig. |21represents the exchange Coulomb energy and is ignored since the typical 
momentum running through the Coulomb line of the order ft compared that in the first diagram, ttid cteft- It is 
shown in appendix C that the Goldstone fields do not generate electric charge distributions. 
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FIG. 2: a is the Coulomb Energy in Eq. 17711 . and b is the ignored exchange Coulomb energy. 
The inhomogeneous block of the stability matrix element, (|68(l . becomes now 

n (k) ^ ( ] = ( 1 + ^*(^)^(^) (82) 

" ^dH*{k)dH{k)^n, ^dH*{k)dH{k)K. k^+ml{k) 

K*{k)n{k) 



A;2 + m|,(A:) 



The stability of the system with respect to the Higgs field requires that linik) > for all k. 

The form factor K{q) and the momentum dependent Debye mass square can be calculated explicitly at weak coupling, 
i.e. 6^ « p, and k « jl. We find that 

n{k) = ^l-Jl^K{k\5p) (83) 



with 



and 



K{k\5^l) - zAtV f dx- ^ ^, \, , , , (84) 



with the first term the Debye mass of the normal phase and I(q\5fi) the function defined in the section III. 
Zero momentum limit: 



Notice that 

k^ov^oc \dH*{k)dH{k)KQ ~ V9A2/^Q 



lim lim ( = (^) , (86) 



^1(^)-^HP}. (87) 



and 



We have 

k^O 



fimf .(^) , (89) 

- ^ dH*{k)dH{k)^ V9A2y^, 



and the charge neutrality stabilize also the inhomogeneous Higgs field with the momentum much smaller than the 
inverse coherence length and the inverse Debye length. 
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FIG. 3: ^£r(red dashed-dotted line) and Ah (green solid line) are plotted as functions of A/5fi in the g2SC phase, and the 
black dashed line shows the corresponding quantity in the 2SC phase. 
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FIG. 4: The function II_ff(fe) (red dashed line), n_ff(fc) (black solid line) and the Coulomb energy (green dash-dotted line) as 
functions of scaled-momentum k/A, in the case of Sfi = 2A and (e^^2)/(47rA^) = 1. 



In the static long- wave length limit, we have 



ml{Q) = ' + (90) 
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and 



4e&^^ 



(91) 



The positivity of the Debye mass square is ensured by the positivity of the M- matrix defined in Ht)2|) . It follows from 
Eqs. H39|) . H83|l . (|90|l and (|91|) that the Higgs self-energy including Coulomb energy correction 



Ah = ilniO) = f 



(92) 



is always positive for the whole range of g2SC state. It means that the Sarma instability in the gapless phase can be 
cured by Coulomb energy. This is shown in Fig. |31 where the red dash-dotted line indicates the Ah and the green 
solid line indicates Ah- 

Nonzero momentum case: 



While the Sarma instability in g2SC phase can be cured by Coulomb energy under the constraint of charge neutrality 
condition, it is not sufficient for the system to be stable, even if the chromomagnetic instabilities are removed, say 
by gluon condensation. One has to explore the Higgs instability by calculating the self-energy function Il{k) in the 
whole momentum space, which amounts to value the three basic functions /(fc|(5/i), J(fc|(5/i) and K(k\dfi) defined in 
(I47|l . (|48|) and 1)84(1 . At T = 0, the summation over Matsubara frequencies becomes an integral and we have 



I{k\Sfi) = — J 



J dz 



A>+o+ ^ y/z^ + A^ (z2 + A2 + ifc2a;2 



(93) 



J{k\Sn) 



A2 
2 



— / dxx^ ( 

TT J-1 V 



and 



■i(5^-0+ 

-oo-\-i6fi 



00+2(5 /J, 



A>+0+ ^'^''Vz2 + A2 (z2 + A2 + ifc2a;2^ 



(94) 



dz 



V^2 + A2^z2 



A2 + ifc2a;2 



(95) 



Further reductions of these integrals are straightforward but rather technical. The details are deferred to Appendix 
IbI Here we present only the numerical results in Fig. 2]and Table I. 

Fig.^shows the Higgs self-energy H//(fc) (the red dashed line) the Coulomb corrected Higgs self-energy tlnik) (the 
black solid line) and the Coulomb energy Ecoui (the green dash-dotted line) as functions of scaled- momentum fc/A, 
in the case of dn = 2A and {e^p.'^)/{ATrA'^) = 1. 

Eventhough the Higgs instability can be removed by the Coulomb energy for small momenta, it returns for in- 
termediate momenta. This phenomenon persists for a wide range of gap magnitude, < A < 0.866^/i and for all 

-2 

strength of the Coulomb interaction, measured by the dimensionless ration 77 = . Within the narrow range 
0.866(5/x < A < Sfi, the Higgs instability could be removed if the Coulomb interaction were sufficiently strong. For 
rj = 1, We found that flnik) > for all k if 0.998(5/i < A < In terms of the values of the parameters of NJL 
model, we have aep.^ < S^. Therefore the electric Coulomb energy cannot cure the Higgs instability for a realistic 
two flavor quark matter. 

Negative H^f(fc) indicates the Higgs mode is unstable and will decay [s^. It is noticed that H/f (fc) reaches its 
minimum at a momentum, i.e., k ~ 4 A, which indicates that a stable state may develop around this minimum, we 
characterize this momentum as kmin ■ The inverse fc~ is the typical wavelength for the unstable mode |3^. If mixed 



phase can be formed, the typical size I of the 2SC bubbles should be as great as i.e., I ~ which 

turns out to becomparable to the coherence length of 2SC in accordance with Eq. (|38|) . (In the case of homogeneous 
superconducting phase, kmin = 0, and / 00.) Considerin g th at the coherence length ^ of a superconductor is 
proportianl to the inverse of the gap magnitude, i.e., ^ ~ A~^ [sj, therefore, a rather large ratio of kmin/ A means a 
rather small ratio of l/^. When l/^ < 1, a phase separation state is more favorable [s^- 

In Tabled we show the Coulomb corrected Higgs self-energy at zero momentum n^f(O), its minimum, the corre- 
sponding kmin/ A and fcmin/Ao ( with A and Aq related by eq. H38|) for different parameters of A/(5/i and aeA'^/fi^. 
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TABLE I: The minimum of the self-energy of Higgs field for various mismatch parameters and strengths of Coulomb interaction. 

It is found that the momentum for the Higgs self-energy to reach minimun is not sensitive to the electromagnetic 
coupHng strength, rather very sensitive to the ratio of Sfi/A. Larger 6fi/A is, larger kmin/A will be. But the ratio of 
this momentum to the gap of 2SC under the same pairing strength is close to one for all cases being examined. 

The reader has to keep in mind that our result of the Higgs instability only indicates some kind of inhomogeneous 
states whose typical length scale is comparable to the coherence length of 2SC. Further insight on the structure of the 
inhomogeneity cannot be gained without exploring the higher order terms of the nonlinear realization (|25|l . A more 
direct approach to obtain the favorite structure of the ground state is to compare the free energy of various candidate 
states, which include the mixed phase, the single-plane wave FF state, striped LO state and multi-plane wave states. 
The Coulomb energy and the gradient energy have to be estimated reliably. We leave this analysis as a future project. 

VI. SUMMARY AND OVERLOOK 

To understand the pair-breaking superfluidity or superconductivity induced by mismatched Fermi momenta, one 
has to go beyond the standard framework of BCS theory. In this paper, we offered a nonlinear realization framework 
of the gauged Nambu-Jona-Lasinio model, where the fluctuations of the Higgs and pseudo Nambu-Goldstone fields 
can be investigated simutaneously. 

It is found that, with the increase of mismatch between the Fermi surfaces of the pairing quarks, not only Nambu- 
Goldstone currents can be spontaneously generated, but the Higgs filed becomes unstable in the gapless phase. As 
we already knew, the instability of the Nambu-Goldstone currents naturally induce the formation of the single-plane 
wave FF state. In this paper, we focused on the Higgs instability which has not drawn much attention yet. 

The Higgs instability is the origin of spatial inhomogeneity, and is not prohibited by global constraints. In the 
case of the imbalanced neutral atom systems, the instability extends from arbitrarily low momenta-albeit not zero-to 
momenta much higher than the inverse coherence length. Without a long range force a macroscopic phase separation 
is likely to be the structure of the ground state since it minimizes the gradient energy of the inhomogeneity without 
additional costs. The thermodynamics of the phase separation can be obtained following the Gibbs construction of 

For the case of two flavor quark matter that is globally neutral, the instability is absent for high momenta and can 
be removed by the electric Coulomb energy for low momenta. But it remains in a window of intermediate momenta. In 
the regime where (5/i >> A, the wavelength for the unstable Higgs mode is alway comparable with the gap magnitude 
of 2SC phase, which indicates the tendency to form 2SC bubbles. Because of the Long range Coulomb force, a 
macroscopic phase separation is prohibited while a crystalline structure is more favorable. Possible candidates include 
the heterotic mixture of normal and superconducting phases and the multi-plane-wave LOFF state. 

Although our calculations are limited to the case of g2SC, the existence of the Higgs instability in the absence of long 
range interactions is generic for all gapless superfluidity/superconductivity that are subject to the Sarma instability. 
Our general formulation can be readily applied to the system with several invariant order parameters and constraints. 
It would be interesting to examine whether the electric and color Coulomb energies are capable of eliminating the 
Higgs instability completely in the gCFL phase. 

It is quite surprising that the Higgs instability shows up at some finite momentum, though it is removed by the 
constraint of the charge neutrality. Following the discovery of the chromomagnetic instability of g2SC at T = 0, the 
authors of _9j examined the Meissner mass square at T 7^ and found that the instability goes away for T sufficiently 
close to Tc- In the light of our result of the Higgs instability, this criterion is not sufficient for the stability. One 
has to extend the analysis of |^ to all momenta of the chromomagnetic polarization function as well as that of the 
Higgs self-energy function at T ^ in order to map out the stability domain of g2SC in the phase diagram of the 
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two flavor quark matter. This amounts to value the functions I(k\Sfj,), J{k\Sfi) and K(k\Sfj,) at nonzero T, which is 
straightforward. We hope to report our progress in this direction in near future. 
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APPENDIX A: CALCULATION OF Hnik) AND Unik) 

The calculation of the one- loop self energies of the inhomogeneously fluctuating fileds introduced in Section ITTll and 
13 will be illurstrated in this appendix with that of the Higgs field, given by 



(Al) 



where we exihibite the dependence of the propagator on momentum and Matsubara energy separately. It is convenient 
to tranform the propagator into the representation of Ref. 21] via 

Sm{ui,P) = US{lj,P)U^ (A2) 

where U = p+ + p-T2 with pj- = i(l ± pa). It follows that 

S{uj,p)^^ = iuj^o - 7 • P + M70/33 - SfijoTs + A75e^p2, (A3) 
where we have suppressed the chemical potential associated to the 8th color. Correspondingly, 

^"^^^ 7feTre3p25K,p + fc)e^p25K,p), (A4) 

where the matrix under the trace is diagonal in isospin space and the trace with respect to iospins becomes the sum 
over the eigenvalue of T3. 

With respect to the color indexes r, g and b, 

^'-{'7 1) (A5) 



0^ 

with (72 the second Pauli matrix in the color SU(2) subspace. Accordingly, 

^"^^^ = ^-\^T.j (0Tr.V25^^(^)K,P + fc)eV25^^(2)K,p), (A6) 

with ti^g (jqIoj. gu(-2) block of the propagator. The expression of S^^^"^' can be obtained from (|A3|) with the 

replacement of by ia2. The trace in color space extends only to r and g indexes now. 
On writing 

[5S(7(2)^^^^j-l ^ ry^p^[p^{ii^ -d-p) + p.- SPP3T3 + iA7o75Cr2Pl], (A7) 

we notice that the matrix (...) can be projected into two orthogonal subspace via 

A± = i(l±p3a-p) (A8) 
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with a = 707. We find that 

(zcj - a • p)pz + p.- 6HP3T3 + i A7o75Cr2Pi 
= A+(iwp3 -p + Ji~ Sp.psTs + iA7o75(T2Pi) + A_(iwp3 + p + fi - S^ip^r^ + zA7o75(72Pi). (A9) 

Under the weak coupling approximation, uj << p, \p — /i| << p, A << /I and Sp << /2, The contribution of the 
second term to the inverse can be ignored. We obtains that 

<,S!7(2). -A \ r- , - X 1 -A ^-l 70P3 -l-p{iuJ- SpT3)p3 +p-fl + A7o75g2Pl 

iu,p) ^ A^i^u;p3-p+p-Spp3r3+^AJoToa,p,) = + ^5 pr^y + {p - p^ + ' ^^^^^ 

Substituting Eq. IjAlOp into Eq. IIA6|l and carried out the trace over Dirac indexes, NG indexes and SU(2) color 
indexes, we end up with 

(.^ = J_^.r i^n + ^SpT3)' + e^'A^ 

' 2Gd ^^^J {27T)^{u;n + tdpT3)^+^'^ + A^][{u;.^+i5pT3y+e + A^] 

^J__.rS^fJ^ {u^n + iSp)' + C(-A^ 

2Gd {2nf[{iu„ + iSp)^+^'^ + A^][{Lj„+iSp)^+e + A^] ^ ' 

where £, = P — p and ^' = p' — p with jf — p + k. The 2nd equality follows from relabeling the index of the sum 
over n such that w„ ^ — w„. In deriving IjAlip we employed the approximation p ~ p' . The integral over p can be 
approximated according to 

//" 1 />oo 
d^p~2TTp d cose I (A12) 
J-1 J -00 

and the ^-integration can be carried out with residues ( notice that i^' — ^ ~ /c cos 9 ) . We have 

^ n.^^ 1 2/i^r^ 1 ^ [uJn + iSpY 

''"^''>-2Gd tt 4r + ^SpP + A2 '^'^ (ivn + i5pY + A2 + 

where we have denote x = cosO for the integration variable. Substituting in the gap equation 



' ^^'^E . ' -''^^T. , ' I - (A14) 



2rr 1 /•! 



we obtain that 

and this derives Eq. (|46|l . 

APPENDIX B: CALCULATION OF INTEGRALS 

In this appendix, we shall calculate the functions I(k\Sp), J{k\6p) and K{k\Sp) at T — and for Sp > A We have 



dxf / + / Idz ^ (Bl) 



J{k\Sp) = — / dxx^i + ]dz ^ (B2) 

2n J_i ^J-oo+^Sf, JiSi,+0+ V^2 _^ A2(z2 + A^ + \k'^xA 
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FIG. 5: The contour for ^-integral. 



and 



K{k\6fi)^— dx( + )dz ^. (B3) 

Because of the similarity of these integrals, we shall focus on I{k\6ii) below. 

Notice that the ^-integral is interupted by the branch cut from iA to infinity along the imaginary axis ( the other 
runs from —iA to infinity ) and ReV + A2 > throughout the cut z-planc. 

Alternatively, the z-integral can be regarded an integral along a contour running from —oo + i6 to oo + iS and going 
around the branch point iA counterclockwisely plus an integral from iS/i + 0"*" to i6fi — 0+ going around the branch 
point clockwisely. The contour of the former can be shifted to the real axis, leaving the result independent of (5/z. We 
obtain that 



with 



I{k\5n)=I{k\0)+Ii{k\5fi) (B4) 



A2 r°° 1 

^(fc|0) = 7^ / dx dz - , ^ (B5) 

27r V^2 + A2fz2 + A2 + iA;2x2 



and 



a2 [\ piSp.-0+ 1 

Ii{k\Sfi)^— dx dz ^. (B6) 

J_i /i5p+o+ Vz2 + A2(z2 + A2 + ifc2a;2] 



Using the formula 
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and making the transformation for the x-integral, x = ^ sinhS, we end up with 



J(fc|0) 



de- 



4A 

Trfc sinh 
8A 
Trfc 



8A \TT^ 1 



2^^,hltanh^ - Li2(e-«0 + ^1^12(6-2"=) 



and 



,,,,,, 2A /■^'^ ^6* , |V'^^l2-rA2cosh6'- (5^sinh( 
/i(fc|d^) = — / . , „ hi ■ " 



Trfc Jo sinh 6* i/(5^2 _ ^pg]^ 6* + (5^ sinh 9 ' 



where is defined by sinh^?c = Li2(z) is Spence function, i.e. 



(B8) 



(B9) 



Li2(z) = dt 



t 



,2 • 



(BIO) 



The integration of J{k\5^) can be carried foUowing the same procedure and the result is elementary, 
16A3 



J{k\5^i) 



+ 



Trfc-^ 
8A3 



, coshf?r — sinh 6* r 



Trg-^ 



, , , I JSa^ - A2 cosh - Sa sinh 9^ I Su , iJSn^ - - A sinh 6*^ 
cosh 9r In In ■ 



v/(5^2 _ cosh 9c + Sn sinh 0^ | ^ ^^'5/^2 _ ^2 + A sinh 61^ 



(Bll) 



Upon decomposing K{k\6iJ.) into i^(fc|0) + Ki{k\6fi) in analog to /(fc|(5/i) we observe that K{k\0) = and the 
integration of Ki{k\6fi) yields 



4A 

Ki{k\Sn) ^ — Rc 
Trfc 



Li. 



A'?,A„2_ A2y 4 'V 



(B12) 



The small fc expansion of /(fc|(5/i), J{k\Sfi) and -R'(fc|(5/i) can be obtainede easily from their integral representations 
as is shown in the text. The large k behavior of /(fc|0) and J{k\Sfj,) cab be extracted readily from their explicit form 
with >> 1. As to /i(fc|(5/i), we write 



^6*0 = / ^6*0 



d9{). 



(B13) 



The large k expansion can be carried out with the aid of the formula ( the derivation is outlined at the end of this 
appendix ). 



d9 J bp? - A2 cosh 6* + 5^ sinh 61 Tr^ 
■In ■ — 



sinh 6* I ^,5^2 _ ^2 cosh 9 - Sp sinh 61 1 2 



(B14) 



and the expansion of the integrand of the 2nd integral on RHS of (|B13|) according to the powers of e . The large fc 
expansion of K{k\6fi) follows from the formula 



for z > 1. We obtain that 



ReLi2(z) = ^ - iln^ z - Li2f- 
3 2 \z 



(B15) 



r(fcM.-L(- 



21n^-2 + ln^^i±V^^ 
A Sp - W - A2 



(B16) 



J(fcM.-±(: 



4 /21n^_2 + ln^^i±V^!H^ 



A2 



(B17) 
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and 



K{k\6^l) ~ (B18) 
k 



Finally, let us sketch how the formula l)B13|l can be derived. Introduce a function of z, 



de , l + ztanh6l ,^ , 

/W= -Tiling— -TTI (B19) 







sinh 1 — z tanh ^ 



the LHS of (|ljl3() equals the real part of its value at z = — ^=M==. For z < 1, the integrand can be expanded 

according to the powers of z and the integration of the coefficients can be carried out explicitly. The resultant power 
series can be recognized to be that of 7rsin~^ z. The analytic continuation to z > 1 gives rise to the RHS of (|ljl3|l . 

APPENDIX C: ELECTRIC CHARGE DENSITY INDUCED BY GOLDSTONE FIELDS 

On writing the electric charge density induced by Goldstone fields 

Sp{k)=Y,>^a{k)Mk), (CI) 

a 

we have 

Kaik) = --eT^ J j^--^Tl-JoQSM{i^7i,P+ k)j5GaSM{LLln,P), (C2) 

where Q = p3{a + bT3), the matrix G4 = e^pi, G5 — e^p2, Gq — e^pi, Gy — e^p2 and Gs — e^pi for the five Goldstone 
fields. 

With respect to r, g, and b color indexes, the propagator assumes the block structure 



while the matrix {p — 1,2) takes the form 

/ n 1/ ^9 Y n \ 

(C4) 



xp(2 X 1) 
-Xt(l X 2) 

Therefore the trace with respect to the group indexes vanishes for and e^. Consequently 

Ki{k) = K5{k) = Keik) = K7{k) = 0. (C5) 

Coming to K8(fc), we firstly take the transpose of the matrix under the trace on RHS of ljC2p . then we employ the 
identity 

7oGiS^(po,p)G'''7o = -5a/(po, -p), (C6) 
and transform the integration variable from p to — p — fc. We find 

K8(fc) = -«8(A:)=0. (C7) 
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